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Summary 

Elementary  initial-boundary-value  problems  are 
solved  for  (i)  a  plane  flux  compression  generator  with 
field  losses  by  diffusion,  (ii)  a  cylindrical  flux  com¬ 
pressor  with  electromagnetic  radiation,  and  (iii)  elec¬ 
tromagnetic  pulse  generation  by  means  of  magnetic  flux 
pushers.  Due  to  the  moving  boundary  conditions,  the 
analytical  treatment  of  these  problems  leads  to  inte¬ 
gral  equations  and  infinite  systems  of  coupled  differ¬ 
ential  equations,  respectively. 

Introduction 


Initial-boundary-value  problems  (IVBP)  for  the 
compression  of  the  flux  of  typical  magnetic  field  con¬ 
figurations  by  accelerated  liners  are  treated  by  space- 
and  time-dependent  eigenfunctions,  which  take  into  ac¬ 
count  the  temporal  change  of  the  field  space  due  to  the 
motion  of  the  conductor  interfaces.  The  derived  solu¬ 
tions  provide  physical  insight  into  the  phenomena  of 
flux  compression  and  dilution,  flux  diffusion  into  ac¬ 
celerated  conductors,  and  the  occurance  of  electromag¬ 
netic  (EM)  radiation! >2  during  flux  compression.  Fur¬ 
thermore,  the  generation  of  EM  wave  pulses  by  means  of 
accelerated  flux  pushers  is  discussed. 


Fig. 1:  EM  induction  model  with  accelerated  conductor  at 
s(t)  s(t)+a,  and  fixed  liners  at  x=-b  and  x=+d. 

Expansion  of  B(x,t) /B0 =  8(£ ,x)  in  space- and  time- 
dependent  eigenfunctions  of 

£=[x-s(t)/a],  x  =  Kt/a2  ;  0  1  ,  0  ^  x  x  ,  (9) 

gives  the  solution  of  Eqs.  (1)  -  (4)  in  dimensionless 
form 

CO  T 

B(£,x)  =  1  +  2tt  I  n/exp  [-n2r2  (x  -  x')  ]  {  [Bp  (x')  -  1] 
n=l  0 


Plane  Flux  Compression  and  Dilution 
with  Diffusion  Losses 


-  (— 1)” [ B2 (t'  )  -  l]}dx'  •  sin  mrS  .  (10) 


Figure  1  depicts  an  EM  induction  system  with  ide¬ 
ally  conducting  (a  =  =»)  walls  at  x  = -b  and  x  =  +d,  in 
which  a  slab  of  finite  conductivity  (a  <  ■»)  occupying 
the  space  s(t)  <  x  <  s(t)+a  is  accelerated  across  a 
homogeneous,  external  magnetic  field  Bo  =  {0,Bo,0}  to  a 
velocity  v  =  (s(t)  ,0,0}  ,  where  s(t=0)  =  so,  -b<s0<+d. 
The  motion  s(t)  %  0  £f  the  piston  induces  (i)  EM  fields 
E = {0,0,E(x,t)}  and  B = {0,B(x,t) ,0)  in  the  conductor  at 
s(t)  <  x  <  s(t)+a,  and  (ii)  EM  fields  in  the  spaces  be¬ 
hind  [Ep (x, t) ,Bp  (t) ]  and  in  front  [E2(x,t)  ,B2(t)  ]  <jf  it 
(Fig.l).  If  the  electric  displacement  current,  3tE/3t 
is  neglected  as  usual,!  the  magnetic  field  B(x,t)  in 
the  accelerated  conductor  is  determined  by  the  parabo¬ 
lic  IBVP  with  moving  boundaries  at  x  =  s(t)  and  x=  s(t)+a: 


B(x,t=0)  =  B0 

B  (x=s  (t)  ,t)  =  Bp  (t) 

B(x=s(t)+a,t)  =  B2  (t) 
where 

d{  [s(t) +b]Bp(t)}/dt  = 
Bp (t=0)  =  B0 
d{  [s  (t)  +  a-  d]B2(t)}/< 
Bo(t=0)  =  B 


3  x2  ,  s(t)<x<s(t)+a  , 

(1) 

,  s0  <  x  <  s0+a  , 

(2) 

,  0  <  t  <- 1 

(3) 

,  0  <  t  <  t 

(4) 

k9B (x=s (t) , t) /3x  , 

(5) 

> 

(6) 

=  ic3B(x=s(t)+a,t)/3x  , 

(7) 

> 

(8) 

K  =  l/ya,  and  t  is  the  flux  compression  _time.  Equations 
(5)  and  (7)  result  from  the  continuity  nx[E]=l5  of  the 
tangential  electric  field  at  the  rear  (1)  and  front  (2) 
interfaces  of  the  conductor,  where  the  conductor  fields 
have_been  eliminated  by  means  of  Ohm's  law,  E  =  -v*B 
+  kV*B,  and  the  EM  fields  in  the  gaps  (1,2)  have  been 
determined  from  V*E  =  -3B/3t  and  VxB  =  6. 
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Substitution  of  Eq.  (10)  into  Eqs.  (5)  and  (7)  de¬ 
monstrates  that  the  gap  fields  Bp(x)  =  Bp(t)/B0  andB2(x) 
=  B2(t)/B0are  determined  by  the  coupled  Volterra  inte- 
gro-differential  equations  [a(x)  =  s(t)/a,  b  =  b/a,  d  =  d/a] 


[a(x)  +6]Bp(x))  =  B2(x)  -  Bp(x)  (11) 

T  T 

-^r/K(x-x')[Bp(x')-l]dx’  +  ,-f/ L(x-x')[B2(t’)  -l]dx'  , 
UT  0  0 

4t{[cf(t)  +  1-c<]B2(x)}-B2(x)  -Bp(T)  (12) 

T  /  X 

+  ^/k(t-x’)[B,(x')  -l]dx'  r  /l(x  -t'  )[Bp(x')  -1]  dx'  , 

QT  0  l  0 

with 

Bp(x=0)  =  1  ,  B2  (x=0)  =  1  ,  (13) 

where 

00 

K(x-x')  =2  £  exp  [-n2ir2  (x  -  x' )]  , 

n=1  (14) 

CO 

L(x-x')=  2  £  (-l)nexp  [-n2it2  (x  -  x'  )  ]  , 

n=l 


are  singular  kernels  of  convolution  type.  In  view  of 
the  extreme  behavior  of  K(x  -  x' )  and  L(x  -  x' )  for  x'  =  x  , 
the  main  contributions  to  the  integrals  arise  for  t'  St. 
For  these  reasons,  Eqs.  (11)  -  (13)  have  the  approximate 
analytical  solutions : 

x  x  x’ 


Bp(x)  -1- 


exp 


[  -JqH)  (x')dx'  ]  jna  (x' )  exp  [+/n<D  (x  )  dx  ]  dx' 


-'0 


0 


o(x)  +  b  +  p(x)  -  q(x)hp(x)/h2(x) 


B2(x)  =  1  +- 

where 
m(x) 


exp[-f  ui(x')dx'  ]/  d(x')exp[+/iu(x)dx)dx' 

IQ. _ 0,  .  .  .  ..  ..  0 - 

d-  1  -  a(x)  +  p  (x)  -  q(x)h2(x)/hp(x) 
b  +  d-  l  +  2rp(x)  +  q(x)  ] 


hp(x)h2(x)  -  q(x){b  +  d-  l  +  2[p(x)  +q(x)]} 


(15) 


(16) 


,  (17) 
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hi  (x  )  =  a  (t  )  +  b  +  p  (x  )  +  q  (t  ) 
h2(x)  =  d  -  1  -a(x)  +p(x)  +  q(x) 


(18) 


r 


and  ^  m 

P(t)  =  +/k(t  -  t'  )dx'  [1  -  4  I  4^P(-n2lr2T)]  >0  ,  (19) 
0  3  71  n=ln 

t  °°  n 

q  (x  )  =  -/l(t  -  x')dt'=i  [1+  I  ifaphV,)]  >0  , 

0  n=1  (20) 

The  results  (10)  and  (15)  -  (16)  are  valid  for  diffusion 
times  tD  -  yaa2  >>  a/c  (since  the  electric  displacement 
current  has  been  neglected)!  an(j  nonrelativistic,  but 
otherwise  arbitrary  piston  speeds  s(t)  =  ad(T). 

As  distinguished  from  ideal  flux  compression 
theory,  B^(x)  and  B2(t)  remain  finite  for  complete  im¬ 
plosion  t  +  t  of  the  left  gap  [o(x) -»■  -  fa]  or  right  gap 
[a  (x)  d  -  1] ,  respectively  (Fig.  1),  since  p(x)  >  0  and 
q(x)>0.  Furthermore,  B^(x)#l  and  B2 (x)  %  1  for  d(x)  %  0 
[compression  (>)  and  dilution  (<)  of  flux],  and  %t2(T) 
=  1  for  d(x)  =0.  The  current  induced  per  unit  width 
Ay  =  1 m in  the  moving  conductor  is 

I(t)  »  (B0/p)I(x)  ,  I(x)  =B2(t)  -B1(x)  ,  0<x<?  .  (21) 

As  an  illustration,  consider  a  flux  compression 
generator  with  CT=6xlO^S/m  (Cu)  ,  a=10“2m,  b  =  d~10_-*-m, 
sG  =  0  m,  s  (t)  =  v0H(t)  ,  vo  =  103m/sec,  t  =  (d  -  a) /v0  ~  10“^ 
sec.  In  this  case,  Eqs.  (15)  and  (16)  give  Bi(t)  -  B0/2 
and  B2(t)  ~  1()3bo.  The  observed  maximum  field  will  be 
smaller  due  to  nonlinear  diffusion  losses,  which  become 
important  for  intensities  B2  i  Bcr =  (ycvT)3/3  =  6. 5x10  T 
for  Cu  at  T  =  103“K. 

The  order-of-magnitude  of  B2(t)  and  B2(t)  does 
not  change  if  the  end  planes  x  =  -b  and  x  =  +d  are  formed 
by  thick  slabs  of  the  same  conductivity  a < “.  For  ap¬ 
plications,  it  is  noted  that  the  a=“  planes  x=-b  and 
x  =  +d  can  also  be  interpreted  as  symmetry  planes  for 
systems  with  four  moving  liners  of  conductivity  o. 

Axial  Flux  Compression  with  EM  Displacement 

In  flux  compression  experiments  EM  wave  fields 
are  observed, ^  which  are  due  to  the  displacement  current 
3eE/3t.  The  implosion  of  a  cylindrical  liner  (of  con¬ 
ductivity  a  =  °°  and  initially  at  R(t=0)  =  R0,  Fig.  2)  with 
a  speed  R(t)<0  against  an  initially  homogeneous  magne¬ 
tic  field  B(r,t=0)  = {0,0,Bo}  excites  by  electromagnetic 
induction  the  wave  fields 

B(r,t)  =44  [rA(r,t)]ez  ,  E(r ,t)  =  -  eQ  *  (22^ 

The  azimuthal  vector  potential  A(r,t)  is  determined  by 
the  hyperbolic  IBVP  with  moving  boundary  at  r  =  R(t): 


32A_  2rl  3  ,JAS  A 


3t2 


r3r  '  3  r'  r2 
1  - 


A(r,t=0)  =  ^  V 
3A(r  ,t=0)/3t  =  0 
A(r=0,t)  =  0 
A(r=R(t),t)  =i-B0R2/R(t)  ,  0<t<t 


R  -  i 
0  % 


0  <  r  <  R(t) 

9 

(23) 

0  <  r  <  Rq 

9 

(24) 

0  <  r  <  Rq 

9 

(25) 

0  <  t  <  t 

9 

(26) 

0  <  t  <  t 

9 

(27) 

le  hypothetical 

time 

t  of 

complete  implosion,  R(t=t)  =  0.  Equations  (24)  ,  (26) 
and  (27)  specify  that  B(r,t=0)  =  Bo,  E(r,t=0)  =  0,  and 
R(t) 

J  B(r,t)2irrdr  =  itR2B0. 

0 


r 


Fig.  2:  pM  fields  E(r,t)  and  B(r,t)  induced  by  imploding 
liner  at  r  =  R(t). 


Equation  (26)  holds  since  the  flux  within  the  circle 
r = 0  is  necessarily  zero.  The  hyperbolic  wave  fields 
i|;(r,t)  are  separated  from  the  quasi-static  elliptic 
compression  fields  by  means  of  the  ansatz 

A(r,t)  =  iHr.t)  +-|B0[R0/R(t)  ]2r  ,  0<r<R(t)  .  (28) 

Equation  (28)  transforms  the  IBVP  (23)  -  (27)  into  an 
IBVP  with  a  source  field  Q(r,t)  =  B0[R'(t)/R(t)  -  3R(t)2/ 
R(t)2] [R0/R(t) ]2r  and  homogeneous  boundary  conditions, 
iKr=0,t)  =  0  and  iJ)(r=R(t)  ,t)  =  0.  For  the  latter  reason, 
ijj(r,t)  tan  be  expanded  in  the  r-  and  t-dependent  eigen¬ 
functions  fk(r,t), 

OO 

i])(r,t)  =  2[R(t)]"2  ^  'IV  (t)  [J0  (<%>  ]-2f]j(t ,  t)  , 

R<t,  “ 

'l'k(t)  =  /  i}i(r,t)fk(r,t)rdr 
where  ^ 

fk(r,t)  =  J1(alcr/R(t))  ,  J1(ak)  =  0,  k  =  l,2,3,,,,“  . 

(30) 

The  eigenvalues  ak  follow  from  the  boundary  condition 
tKr=R(t)  ,t)  =  0,  whereas  iKr=0,t)=0,  since  J2(0)=0. 

By  means  of  this  extended  Fourier  method,  the  solution 
for  i|j(r,t)  is  found  which  gives,  under  consideration  of 
Eqs.  (22)  and  (28),  the  EM  field  solutions  in  the  im¬ 
ploding  cylindrical  cavity  0<r<R(t): 


*k(t) 


“kr 


A(r ’ t)  2B°[R(t)]  r  +  R(ty2k^1J0(ak)2  Jl(R(t) 


and 


-  Vk(t) 


V 


B(r,t)  MR(t)]+R(t)3J^Jo(ak)2  J°  (R(t) 


)  ,  (31) 


(32) 


E(r  t]  -  R(t)B  [  R°  ]2  r  1 2B(t)  y  v  J 

E(r,t;  R(t)B0  J  k(-t-)  A.  t_  i'™. -12  no'i 


(t)  J  R(t)  TR(t)3ki1Jc(ak)2  R(t)  0  '‘R(t)  ' 


R(t)3 


”  [R(t)ip,  (t)  -  R(t)<k(t)]  a,  r 
l  TTivTTJ  Jl  ^R(t)  ^  *  ^33^ 


k=I  J0(ak)2 


The  Fourier  amplitudes  <[Jk(t)  of  the  hyperbolic  wave 
fields  are  determined  by  the  infinite  set  of  coupled 
differential  equations  [derived  from  the  inhomogeneous 
wave  equation  for  i)i(r,t)]. 


d2t 


k  _  „  R(t)  d^k 


(a,  c) 


dt2  R(t)  dt 


4.  r_L__  1  9'R(t)z, 

[  R(t)2  2R(t)  +  (2  3  ak  R(t)2  ^k 


00  di// . 

’  Qk(t)  -  I  [Mfcj  (t)ipj  +Nkj  (t)  ] 

jrk 


(34) 


with 


110 


ipk(t=0)  =  0  ,  d<|>k(t=0)/dt  =  0  ,  k  =  1,2,3, 


,  (35) 


where 

Qk(t)  =  B0R2{3R(t)2/R(t)  -  R(t)}  J0(ak)/ak 


and 


R(t)2 ,  R(t)  R(t) 
R(t)2  R(t)  R(t) 


4a? 

_ _ _ 

2  2-12-2 
“k-oj  “k  “j 


a  .a  J  (a,  ) 


J„  (a.) 


rw 


•  v  a  .a, 
xt  /  R(t)  j  k 

kj(  }  4  R(t)a2  -a2  1  Jq  (aj  ) 


5  kj  ■ 


,  (36) 

V  • 

(37) 

(38) 


The  coupling  of  the  wave  modes  k  in  Eq.  (34)  is  due  to 
the  boundary  motion  R(t)  ,  since  Mk1  (t)  4  0  and  Ni^  (t)  4  0 
for  R(t)  4  0. 

Equations  (31)  -  (33)  indicate  that  the  complete 
hyperbolic  solutions  are  the  sum  of  (i)  quasi-static  (S) 
fields  and  (ii)  transient  wave  (R)  fields  (Fourier 
series).  By  Eq.  (34),  the  "radiation"  fields  are  rela- 
tivistically  small  since 


'l'k(t)  0  .  R(t)2 / c2  ->■  0  .  (39) 

Accordingly,  flux  compression  theory  without  displace¬ 
ment  current  is  approximately  valid  for  nonrelativistic 
liner  speeds,  R(t)2<<c2.  Since_^k(t)  -  B0  (R(t) /c)2R2x 
R(t)  by  Eq.  (34)  and  B2/2p  =  (Bg -j^Bk^2/2p,  the  radiative 
energy  density  is  of  the  order  Bg-BR/p~  (B2/p)  (R(t)/c)2. 
For  these  reasons,  an  explanation  of  magnetic  flux  com¬ 
pression  during  gravitational  collapse  of  magnetic  stars 
withspeeds  |R(t)js  c  must  be  based  on  Maxwell's  equat¬ 
ions  with  displacement  current. 


EM  Wave  Pulses  from  Flux  Pushers 

Theoretical  investigations  indicate  that  EM  waves 
are  excited  at  the  surface  of  conductors  which  are  ac¬ 
celerated  in  magnetic  fields. 1,2  in  cylindrical  flux 
compression  experiments,  the  emission  of  EM  pulses  from 
the  openings  of  the  cylinder  ends  has  been  observed.^ 
These  explosion  generated  EM  pulses  have  a  low  charac¬ 
teristic  frequency  (Fourier  spectrum)  u)  ~  At-1  since  the 
duration  of  the  explosive  conductor  motion  is  of  the 
order  At  -  lO'^sec. 

As  an  illustrative  model  for  the  inductive  gener¬ 
ation  of  EM  pulses,  consider  a  plane  copper  piston  of 
thickness  d  with  its  front  surface  initially  a £  x=  0, 
which  borders  on  a  homogeneous  magnetic  field  B0  =  {0,Bo,0} 
filling  the  half  space  x>0  (Fig.  3).  At  time  t  =  0, 
this  piston  is  pushed  with  a  velocity  v=  (a(t)  ,0,0)  into 
the  B0-field  region  x>a(t)  >  0.  The  duration  At  of  the 
impulsive  piston  motion  is  assumed  to  be  small  compared 
with  the  EM  diffusion  time  tp  =  p0cd2  -  10-2sec  (for  a  = 
6xl()ls/m  and  d =  10_2m) ,  so  that  the  piston  can  be  trea¬ 
ted  as  a  perfect  conductor  (a  =  °°)  .  For  a  piston  of 
area  AA=lm2,  the  EM  pressure  force  F,  its  work  AW  along 
a  piston  path  Ax=lm,  and  the  generated  power  P  for  a 
shock  wave  driven  piston  with  speed  v  ~  10^m/sec  are 

F  =  AAB^/y  a  106N,  AW=FAx=106J,  P  a  KjiOj/sec 

11  B0  =  (4tt/10)^2  =  i.i  t.  The  duration  of  the  piston 
motion  is  At  -  Ax/v  -  10~^sec.  The  expanded  work  AW  must 
appear  in  the  form  of  an  EM  energy  pulse. 


Fig.  3:  Plane  flux  pusher  at  x=a(t)  with  transverse 
magnetic  field  B  and  EM  pulse  S. 


For  the  plane  magnetig  flux  pusher  in  Fig.  3,  the 
EM  fields  are  of  the  form,  B = { 0,-3A(x, t) /3x,0)  and  E  = 
{ 0,0,-3A(x, t) /3t) .  The  vector  potential  is  determined 
by  the  IBVP  with  a  moving  boundary  condition  at  the 
piston  front  face  x=a(t): 


32A/3t2  =  c232A/3x2  , 

A(x,t=0) = -B0x  , 

3A(x,t=0)/3t  =  0 


a(t)  <  x  <  “ 
0  <  x  <  “ 

0  <  x  <  " 


,  (40) 

,  (41) 

,  (42) 


[3A/x,  t) /3t  +  a(t)  3A(x,  t) /3x]  _  .  ,  =0 

x  att; 


t  >  0 


(43) 


Equations  (^1)  - ^43^  consider  that  B(x,t=0)  =  Bo, 
E(x,t=0)  =  0,  and  E  +  vXB  =  0  at  the  piston  (a=“)  surface 
x=a(t).  The  solution  of  Eqs.  (40)  -  (43)  is  by  the  La¬ 
place  transform  method 


a+i“ 


A(x,t)  =-B0x+^~/_^F(s)es(t  x/c) 


ds 


(44) 


where  the  amplitude  function  F(s)  satisfies  by  Eq.  (43) 
the  integral  equation 


a+i“ 


2tt 


—  /  sF(s)es[t_a(t)/c]ds  =  a(t)[l-a(t)/c]-1B0  . 

(45) 


Equation  (45)  can  be  solved  analytically  for  certain 
piston  motions  a(t),  e.g.,  for  the  Heaviside  step  veloc¬ 
ity  model. 


a(t)=v0H(t)  ,  F(s)  =v0B0(l-v0/c)  ls  2  .  (46) 

From  Eqs.  (44)  and  (46)  one  obtains  th^  s^Lytions  for 
the  EM  fields  and  the  Poynting  vector  S  =  ExB/p: 


v 


A  (x ,  t )  =  -B0x  +  v0B0  (1  -  -fy ~x  (t  -  |)H  (t  - 1) 

,  x  >  v0t  , 

(47) 

E(x,t)=-v0B0(l-^-)-1H(t-f) 

,  x>  v0t  , 

(48) 

B(x,t)=Bo+^(l-^)-lB0H(t-f) 

,  x>  v0t  , 

=  0 

.  x<  v0t  , 

(49) 

and 

S  =  -EB/p  =  S0  +  S 

(50) 

where  2 

B„  vD  , 

S0(x,t)  =  v0—  (1-  — )  ^(t--) 

.  x>  v0t  , 

=  0 

,  x  <  v0t  , 

(51) 

B(x,  t)  =~  vQ  —  (1  -  ~)-2H(t  -  ~) 

,  x  >  vQt  , 

=  0 

,  x  <  vct  . 

(52) 

It  is  seen  that  the  energy  flow  density  S  consists  of  a 
"mixed"  (static-transient)  term  S0  =  -EB0/p  and  a  true 
wave  term  S  =  -EB/p ,  where  S~  (v0/c)S0  if  vQ  <<  c.  The 
hypothetical  limits,  S0->-”  and  S-*“  for  v0^-c,are  re¬ 
markable.  For  the  above  B,-,,  v  =v,  and  At  values,  Eqs. 
(51)  and  (52)  give  (AA=lm2) 

SoAAAt-106J,  and  SAAAt  =  (l/3)xl02J 

since  l-v0/c  =  l.  Accordingly,  the  main  energy  flow  is 
in  SQ,  but  the  true  wave  contribution  S  represents  also 
a  significant,  measurable  effect.  The  efficiency  of  the 
production  of  the  wave  pulse  S  is  of  the  order  v0/c  - 
10“^/3,  which  is  comparable  to  the  efficiency  ~10~4  0f 
the  radiation  production  in  pulsed  x-ray  tubes.® 

By  the  same  method,  the  EM  wave  pulses  generated 
by  means  of  more  complicated  flux  pushers  (Figs.  4,5)  can 
be  evaluated.® 


Ill 


Fig.  4:  Cylindrical  flux  pushej:  at  r  =  R(t)  wij:h  azi¬ 
muthal  magnetic  field  B  and  EM  pulse  S. 

The  radial  EM  pulse  generator  in  Fig.  4  consists 
of  a  cylindrical  liner  of  initial  (external)  ra¬ 

dius  R(t=0)=Ro,  which  is  exploded  to  a  radius  R(t)  >  Rq 
at  time  t >  0.  The  liner  carries  a  current  I (t) =  2rR(t) 
xj*(t)  on  its  surface  r  =  R(t)  with  l(t=0)  ~  10^A.  The 
pulsed  liner  motion  R(t)  transverse  to  the  magnetic 
field  Bg(r,t)  of  I(t)  excites  an  EM  energy  pulse  in  ra¬ 
dial  direction, 

Sr= -Ez(r,t)B9(r,t)/u  >  0  ,  R(t)  <  r  <  Ro  +  ct  ,  (53) 

since  the  induced  electric  field  is  Ez(r,t)  *  -R(t)x 
Bg(r,t)<  0.  The  liner  functions  for  a  limited  pulse 
time  6t  (depending  on  its  initial  thickness  &R0) ,  since 
the  conductor  becomes  too  thin  and  too  transparent  at  a 
critical  expansion  radius  Rcr  =  R(t=At). 


Fig.  5:  Cylindrical  flux  pusher  at_r  =  R(z,t)  with  local¬ 
ly  deformed  magnetic  field  B  and  EM  pulses  S. 


The  axial  EM  pulse  generator  in  Fig.  5  consists  of 
a  cylindrical  liner  (a  ->  °°)  of  initial  (inner)  radius 
R(z,t=0)=Ro  and  contains  an  initially  homogeneous  mag¬ 
netic  field  Bz (r ,z ,t=0)  =  B0  in  the  axial  direction, 
where  B0  ~  IT.  This  liner  is  imploded  with  the  speed 
3R(z,t)/3t<0  within  the  axial  region  -z0<z<+zo.  The 
resulting  compression  of  the  axial  magnetic  flux  induces 
an  electric  wave  pulse 


Eg  (r ,  z ,  t)  =  - 


i  f  3Bz(r»z>t) 


2rr 


3 1 


2irrdr  <  0  ,  0  <  |  z|  <  ct  ,  (54) 


since  3Bz(r ,z,t)/3t >  0.  The  energy  flux  of  the  induced 
EM  fields  is  in  the  axial  directions, 


sz  =  -Eg  (r,z,t)Br(r,z,t)/p  %  0  for  z®  0  ,  0  <  [  z|  <  ct  ,  (55) 

since  Br(r,z,t)%0  for  z%  0.  Thus,  EM  wave  pulses  are 
ejected  from  both  ends  of  the  cylinder  during  magnetic 
flux  compression  4  (Fig.  5). 

The  generation  of  EM  wave  fluxes  S  by  means  of 
flux  pushers  would  be  very  efficient  if  relativistic 
liner  motions  v£  c  were  realizable.  For  this  reason, 
the  transformation  of  chemical  energy  of  explosives,  via 
work  against  the  EM  back  pressure  of  an  initially  (qua¬ 
si-)  static  magnetic  field,  into  EM  radiation  represents 
an  interesting,  but  not  readily  solvable  problem.  If 


Ohmic  losses  are  negligible  (At/tjj-x  0)  ,  this  energy 
transformation  mechanism  follows  directly  from  Poynting's 
theorem,  according  to  which  the  work  3W/3t  of  the  flux 
pusher  per  unit  time  equals  the  EM  energy  flow  S  through 
the  surface  (dli)  of  the  system  (a-*-  «°) , 

3W/3t  =  -u-1$  ExB"  da. 


Appendix:  Derivation  of  Eq.  (34) 

Substitution  of  Eq.  (28)  into  Eq.  (23)  results  in 
the  inhomogeneous  wave  equation  for  i|j  ( r,t), 


f&"c2[r&(r^)~'£1+Q(r’t)  ’  0<r<R(t)  ,  (Al) 
where 


r*R*(t)  „  R(t)2,  r  Ro  1?i 


Q(r’t)=Vf^-3f^][^r 


(A2) 


is  a  known  source  field.  The  Fourier  amplitudes  tk(t) 
are  determined  by  multiplication  of  Eq.  (Al)  with  the 
eigenfunctions  (30)  and  subsequent  integration  over  the 
area  irR(t)2, 

R(t)  ..  ,  .  .  R(t) 

fkrdr=c2/  [i— (r^-)  -i]fkrdr  +  /  Qfkrdr. 

0  0  (A3) 


0  3t2 


Partial  integrations,  under  consideration  of  the  bound¬ 
ary  values  ^(r,t)  =  0  and  fk(r,t)  =  0  for  r  =  0,  R(t),  yield 


?(t)  3  3 

/  t1 

'  r* 

0 


(t|r)  -~2]fkrdr  =  -[“k/E(t)]2'('k(t) 


r  3r  3r 


and 

/  3l2  fkrdr  = 


d2^k(t) 

-/ 


dt2 


?(t)  3jj;_3fk 


32f, 


0 


[23¥TE+,J;3^]rdr 


(A4) 


(A5) 


Equation  (A3)  gives,  under  consideration  of  Eqs.  (A4) 
and  (A5) ,  the  differential  equation  for  ipk(t) : 


d2<p  (t) 

-  ~df ~ ■ -  +  [akc/R(t)]2^k(t)  -  Qk(t)  =  C2k(t)  (A6) 


where 


3  t 


'I'.f.  32  f, 

+(R(t)2)T?_]rdr  • 

(A7) 


The  Bessel  integrals  in  Eq.  (A7)  can  be  evaluated  in 
closed  form.  Thus,  one  finds  Eq.  (34)  from  Eq.  (A6) . 
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